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ABSTRACT 

Few galaxies have been found between the redshift ranges z 1 probed by magnitude-limited surveys 
and z 3 probed by Lyman-break surveys. Comparison of galaxy samples at lower and higher redshift 
suggests that large numbers of stars were born and the Hubble sequence began to take shape at the 
intermediate redshifts 1 < z < 3, but observational challenges have prevented us from observing the 
process in much detail. We present simple and efficient strategies that can be used to find large 
numbers of galaxies throughout this important but unexplored redshift range. All the strategies are 
based on selecting galaxies for spectroscopy on the basis of their colors in ground-based images taken 
through a small number of optical filters: GIZi for redshifts 0.85 < z < 1.15, GIZz for 1 < z < 1.5, and 
U n G1Z for 1.4 < z < 2.1 and 1.9 < z < 2.7. The performance of our strategies is quantified empirically 
through spectroscopy of more than 2000 galaxies at 1 < z < 3.5. We estimate that more than half 
of the UV-luminosity density at 1 < z < 3 is produced by galaxies that satisfy our color-selection 
criteria. Our methodology is described in detail, allowing readers to devise analogous selection criteria 
for other optical filter systems. 

Subject headings: galaxies: formation, galaxies: high-redshift 



1. INTRODUCTION 

As the photons from the microwave background stream 
towards earth they are gradually joined by other photons, 
first by those produced in the occasional recombinations 
of the intergalactic medium, later by those cast off from 
cooling H 2 molecules, later still, after hundreds of mil- 
lions of years, by increasing numbers from quasars and 
galaxies, by bremsstrahlung from the hot gas in galaxy 
groups and clusters, and finally, shortly before impact, 
by photons emitted by the Milky Way's own gas and 
stars and dust. All reach the earth together. They pro- 
vide a record of the history of the universe throughout its 
evolution, but it is a confused record. Two photons that 
simultaneously pierce the same pixel of a detector may 
have been emitted billions of years apart by regions of 
the universe that were in vastly different stages of evolu- 
tion. One of the challenges in observational cosmology is 
to separate the layers of history that we on earth receive 
superposed. 

This paper is concerned with a small part of the chal- 
lenge: locating galaxies at redshifts 1 ^ z ^ 3 among the 
countless objects that speckle the night sky. The red- 

1 Based in part on observations obtained at the W.M. Keck Ob- 
servatory, which is operated jointly by the California Institute of 
Technology, the University of California, and NASA, and was made 
possible by a gift from the W.M. Keck Foundation. 

2 Harvard Society Junior Fellow 

3 Current address: Carnegie Observatories, 813 Santa Barbara 
St, Pasadena, CA, 91101 

4 Current address: Department of Astronomy, University of Cal- 
ifornia, Berkeley, CA 94720 



shift range is interesting for a number of reasons. Dur- 
ing the ~ 4 billion years that elapsed between redshifts 3 
and 1, the comoving density of star formation was near 
its peak, the Hubble sequence of galaxies may have been 
established, and a large fraction of the stars in the uni- 
verse were born (Dickinson et al. 2003, Rudnick et al. 
2003). Observing star- formation at these redshifts will 
likely play an important role in our attempts to under- 
stand galaxy formation. 

Identifying galaxies at 1 < z < 3 would be easy if ob- 
serving time were infinite. One could simply obtain a 
spectrum of every object in a deep image and discard 
those at lower or higher redshift. In practice telescope 
time is limited and one would like to spend as much of 
it as possible studying objects of interest — not finding 
them. Our aim is to present simple strategies that allow 
observers to locate star-forming galaxies at 1 < z < 3 
with only a minimal loss of telescope time to deep imag- 
ing or spectroscopy of objects at the wrong redshifts. 

We shall assume throughout that the goal is to iden- 
tify star-forming galaxies that lie inside a narrow range of 
redshifts within the wider span 1 < z < 3. Concentrat- 
ing on a narrow range has obvious observational advan- 
tages; galaxies at similar redshifts have similar spectral 
features that lie at similar observed wavelengths, and 
consequently the spectrograph can be chosen and con- 
figured in a way that is optimized for every object on a 
multislit mask. It is also necessary for many scientific 
programs. One might want to discover where galaxies lie 
relative to the intergalactic gas that produces absorption 
lines in the spectra of a QSO at z ~ 2.5, for example, 
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or to create a large and homogeneous sample of galax- 
ies at a single epoch in the past that can then be com- 
pared to existing samples of galaxies in the local universe. 
Projects with wildly different goals may benefit less from 
the strategies we advocate. Standard magnitude-limited 
spectroscopy or photometric redshifts might be prefer- 
able if (for example) one wanted to create a sample that 
contained all galaxy types at all redshifts. 

Our approach exploits the fact that the wavelength- 
dependent cross-sections of various common atoms and 
ions give distinctive colors to galaxies at different red- 
shifts, ft has been recognized for decades that mea- 
suring a galaxy's brightness through a range of broad- 
band filters should therefore provide some indication of 
its redshift (e.g., Baum 1962; Koo 1985; Loh & Spillar 
1986). Recent work suggests that a redshift accuracy of 
a z ^0.1 can be achieved for galaxies at £ z < 6 given 
extravagantly precise photometry through 7 filters that 
span the wavelength range 0.3 ^ A 2.2/zm (e.g., Hogg 
et al. 1998; Budavari et al. 2000; Fernandez-Soto et al. 
2001; Rowan-Robinson 2003). An accuracy <r z ^ 0.1 is 
clearly sufficient for finding galaxies at 1 < z < 3, but 
obtaining the necessary deep images through numerous 
filters consumes enormous amounts of telescope time that 
could be more profitably devoted to galaxy spectroscopy. 
We were led to seek specialized color-selection techniques 
that could identify galaxies at 1 < z < 3 even in com- 
paratively noisy images taken through a small number of 
optical filters. 5 

Such techniques can succeed only if they take advan- 
tage of strong and obvious features in galaxies' spectra. 
No feature is stronger than the Lyman break at 9 12 A 
produced by the photo-electric opacity of hydrogen in its 
ground state. Meier (1976) argued that the strength of 
this break would allow high-redshift galaxies to be iden- 
tified in images taken through just three filters, a claim 
that Steidel et al. (1996; 1999) have since confirmed. 
Although the Lyman break itself is not visible from the 
ground at z ^ 3, other weaker features are, and the suc- 
cess of two-color selection at z > 3 inspired us try to 
develop similar two-color optical selection strategies for 
1 < z < 3. Figure shows some of the spectral fea- 
tures that we had to work with. Most obvious, after 
the Lyman-break, is the Balmer-break at 3700A. The 
strength of the Balmer-break can be estimated from the 
model galaxy spectra described in §|21or from the galaxy 
observations described in § EI §S explains how it can be 
used to find galaxies at 1 £ z ,$ 1.5. Franx et al. (2003) 
and Davis et al. (2003) have also used this feature to find 
distant galaxies. At 1.5 H, z ^ 2.5 no strong breaks are 
present in the optical spectra of galaxies, but the lack of 
spectral breaks is itself a distinguishing characteristic of 
galaxies at these redshifts. §§EHH1 and explain. Our re- 
sults are summarized and discussed in § El Together with 
the Lyman-break technique, the selection techniques pre- 
sented here allow the efficient creation of large samples 
of star-forming galaxies throughout the redshift range 
1 £ z < 5. 

5 As readers will notice, the samples created by our techniques 
will hardly differ from those that would be created if standard 
photometric-redshift techniques were applied to the same noisy 
two-color data; our goal is not to disparage photometric-redshift 
techniques but to adapt them to a new regime. 
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Fig. 1. — UnGlZiz colors of galaxies at redshifts < z < 3. 
Shaded curves show the spectrum of a model star-forming galaxy 
(type Im from §[2J at various redshifts. Unshaded curves show the 
transmission curves of the filters used in this paper. 



2. MODEL GALAXY SPECTRA 

Our development of selection strategies began with 
theoretical models of galaxy spectra. At redshifts z ^ 3 a 
galaxy's optical broad-band colors are determined by the 
mixture of stellar types it contains. This is in turn largely 
determined by a galaxy's star-formation history. Galax- 
ies that formed most of their stars recently have spectra 
dominated by bright and hot massive stars, while galax- 
ies that formed most of their stars in the distant past will 
have spectra dominated by fainter, cooler, but longer- 
lived low-mass stars. We considered five model galaxy 
spectra that were intended to span the range of possi- 
ble star-formation histories. Each model spectrum was 
calculated with the code of Bruzual & Chariot (1996, pri- 
vate communication) and assumed that the galaxy's star- 
formation rate as a function of time, S(t), was a decaying 
exponential: S(t) oc exp(— t/r). The adopted values of r 
and assumed time-lapse since the onset of star formation 
for our five models are listed in Tabled Bruzual & Char- 
lot (1993) show that star-formation histories with these 
parameters reproduce the observed spectra of different 
galaxy types in the local universe. Because a wide range 
of star-formation histories can result in nearly identical 
model spectra, we were not concerned that some of our 
adopted model parameters are physically impossible at 
high redshift due to the young age of the universe. Pa- 
rameter combinations that are more plausible can pro- 
duce similar spectra, and in any case our aim was only 
to have model spectra that roughly spanned the range of 
conceivability. Subsequent empirical refinements of our 
selection criteria would compensate for any shortcomings 
in our model spectra. 

We estimated the colors of galaxies at different red- 
shifts by scaling the wavelengths of our template spectra 
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Table 1. Model SED Parameters 



Name 




age a 


E 


1 


13.8 


Sb 


2 


8.0 


Sbc 


4 


10.5 


Sc 


7 


12.3 


Im 


oc 


1.0 



a Gyr 



by 1 + z, applying the appropriate amount of absorp- 
tion due to intergalactic hydrogen (Madau 1995), and 
finally multiplying the result by our filter transmissivi- 
ties. In some cases, mentioned explicitly below, the tem- 
plate galaxies were first reddened by dust that followed 
a Calzetti (1997) attenuation law. For the reddening 
E(B — V) = 0.15 that appears typical for high-redshift 
galaxies (Adelberger & Steidel 2000) the resulting change 
in galaxy color is not large: 0.22 magnitudes from rest- 
frame 1500A to 2000A, 0.16 from 2000A to 2500A, 0.13 
from 2500A to 3000A, 0.21 from 3000A to 4000A. The 
formulae that we used to calculate template galaxies' col- 
ors can be found (e.g.) in § 4.2 of Papovich, Dickinson, 
& Ferguson 2001. 

3. OBSERVATIONS 

Although our initial ideas for color-selection strategies 
were motivated by model galaxy spectra (§ |2Jl, our fi- 
nal selection criteria were determined by the observed 
broadband colors of galaxies at different redshifts. 

The galaxies we used lay within fields observed during 
our survey of Lyman-break galaxies at z ~ 3 (Steidel et 
al. 2003) or within fields chosen with similar criteria. 
U n G1Z images of each field were obtained as described in 
Steidel et al. (2003). In some cases these images were 
supplemented with i band images to increase our wave- 
length coverage. In one field, the HDF-North, a publicly 
released z image from the GOODS survey (Dickinson & 
Giavalisco 2002) gave us photometric coverage of the en- 
tire optical window. Transmission curves for these filters 
are shown in Figure Q] The typical set of U n GlZi images 
was taken in ~ 1" seeing at the Palomar 200-inch tele- 
scope with exposure times (median) of 7.0, 2.3, 1.8, and 
2.0 hours respectively. All reported magnitudes are in 
the AB system. A complete description of our U n GlZi 
observations can be found in Steidel et al. (2003). 

Galaxies were selected for spectroscopic follow-up on 
the basis of their colors. Initially our spectroscopy was 
exploratory, as we sought to establish whether a certain 
combination of colors reliably indicated that a galaxy 
lay within a targeted redshift range. During this phase 
galaxies were observed on multislit masks that were pri- 
marily devoted to the Lyman-break survey. The spectro- 
scopic set-up for these observations (I'M slits, ~ 10A res- 
olution, 4000 < A < 7000A wavelength coverage, ~ 1.5 
hour integration time with LPJS on the Keck I or II tele- 
scopes) is described in Steidel et al. (2003). Later, after 
our initial color-selection criteria had been validated or 
refined, entire masks were filled with objects that sat- 
isfied them, and the spectroscopic set-up was optimized 
to the targeted redshift. For redshifts 1.5 < z < 2.5, 



where our redshift measurements were based primarily 
on absorption lines in the rest-frame far-UV, the dif- 
ference from the Lyman-break survey set-up was slight. 
Most of these spectra were obtained with the blue arm of 
LRIS using a 400 line mm -1 grism blazed at 3400 A and 
1'.'2 slits. The resulting spectra had ~ 5A resolution and 
stretched from 3100 to 8000A. For roughly half of the 
spectra we concurrently obtained spectra with the red 
arm, usually with a 6800 A dichroic and 400 line mm -1 
8500A grating that provided spectra of ~ 5A resolution 
from 6800A to 9500A. At z ~ 1 our redshift measure- 
ments were based primarily on the [OIIJA3727 doublet. 
Here we benefited from redder, higher-resolution grat- 
ings and could tolerate shorter exposure times. 0'.'8 slits, 
~ 3A resolution, wavelength coverage from 6300A to 
8700A, and 1 hour integration times were typical. All 
data were reduced with procedures similar to those de- 
scribed in Steidel et al. (2003). 

Because our goal was to measure the maximum num- 
ber of galaxy redshifts, we would begin observing a new 
slitmask after the previous slitmask's allotted integration 
time had elapsed, even if (as was usually the case) our 
exposures had not been long enough to produce spectra 
that allowed us to measure a redshift for every object. 
In the best conditions we were able to measure redshifts 
for 90% of the objects. In the worst conditions few of 
our spectra were usable. Averaging over all conditions, 
our net spectroscopic success rate was roughly 75%. We 
have been unable to find any evidence that objects with 
identified and unidentified spectra lie at significantly dif- 
ferent redshifts. In the numerous cases where we were 
able to measure a previously unidentified object's red- 
shift by observing it again on a new slitmask, its redshift 
was similar to those of the objects whose spectra were 
identifiable after the first attempt. Nevertheless readers 
should be aware that our sample suffers from some in- 
completeness. Figure El shows our spectroscopic failure 
rate as a function of U n GlZ color for objects whose spec- 
tra were obtained with the blue spectroscopic configura- 
tion described above and in Steidel et al. (2003). There 
is little evidence that the failure rate depends strongly 
on objects' intrinsic colors within the color selection win- 
dows described below. 

4. COLOR SELECTION AT 1.0 < Z < 1.5 

4.1. Balmer-break selection at z ~ 1.0 

Most galaxies have a break in their spectra at A rcs t ~ 
4000A that is produced by a combination of Hydrogen 
Balmer-continuum absorption in the spectra of B, A, and 
F stars and Call H&K absorption in the spectra of F, 
G, and K stars. The relative strengths of the Balmcr 
and 4000A breaks depends upon the mixture of stel- 
lar types in a galaxy — younger galaxies have stronger 
Balmer breaks and older galaxies have stronger 4000A 
breaks — but few galaxies have no break at all (Figure EJ. 

Figure n suggests that this break should give distinc- 
tive broadband colors to galaxies at z ~ 1. At no other 
redshift will a strong break between 1Z and i be accompa- 
nied by a flatter spectrum at bluer wavelengths. Our first 
guess at photometric selection criteria targeting galaxies 
at z ~ 1 was inspired by figure 01 which shows the GIZi 
colors of model galaxies at various redshifts (S|2J) and of 
stars in our own galaxy. Spectroscopic follow-up of ob- 
jects with colors that are characteristic of z ~ 1 galaxies 
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Fig. 2. — The fraction of spectroscopically observed galaxies for 
which we were unable to measure a redshift as a function of U n G1Z 
color. 
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Fig. 3. — The Balmer and 4000A breaks. Younger galaxies, 
represented by the "Im" model spectrum f§[2"l. have a Balmer break 
near 3700A due to Hydrogen Balmer absorption in the spectra of 
B, A, and F stars. Older galaxies, represented by the "Sb" model 
spectrum (§0, have 4000A Call H&K breaks. These breaks give 
distinctive GIZi colors to galaxies at z ~ 1. 



and distinct from the colors of other objects should pro- 
duce a redshift survey consisting primarily of galaxies at 
z ~ 1. These objects lie to the right of the diagonal line 



in figure 21 suggesting 

K-i > QA(G-K) + 0.2, G - K < 1.8 [FN] 

as reasonable criteria for identifying the z ~ 1 galaxies 
in deep GIZi images. These criteria are called "FN" in 
our internal naming convention and may be referred to 
by that name in subsequent publications. 




% - i 

Fig. 4. — Expected locations of stars and galaxies on a GIZi 
two-color diagram. The stars are from Gunn &; Stryker (1983). 
The curved tracks show the colors of model galaxies of different 
spectral types (see § at various redshifts, starting at z = and 
increasing clockwise. Circles mark redshift intervals of dz = 0.1; 
large circles mark z = 0.9, 1.0, 1.1. Dust reddening (Calzetti 1997) 
to E(B — V) = 0.15 was assumed. The arrow indicates the color 
change corresponding to AE(B — V) = 0.2 for an Sc galaxy at 
z = 1.0. One would expect to find galaxies at z ~ 1 and little else 
to the right of the diagonal line. Since the line is nearly parallel to 
the reddening vector, this conclusion does not depend strongly on 
galaxies' assumed dust content. Sb galaxies fall to the right of the 
line even for redshifts significantly larger than 2 = 1, but higher- 
redshift Sbs are not a major component of our sample because Sb 
spectral types are significantly less common than Sc at z ~ 1 (cf. 
figure |FJ and because higher redshift sources are fainter and less 
easily able to satisfy any apparent-magnitude limit. The dotted 
line shows the colors of a model galaxy at z = 1.0 that has been 
forming stars at a constant rate. Crosses mark its colors after 
3,10,30,100,300,1000 Myr of star-formation. The galaxy becomes 
redder as it ages and begins to satisfy our selection criteria for 
t £ 30Myr. 

Even if galaxies at z ~ 1 had colors that matched the 
models perfectly, and even if we suffered no photometric 
errors, figure 0] makes it clear that some galaxies at z ~ 1 
would not satisfy the selection criteria above. These are 
galaxies with extreme star-formation histories. At one 
extreme are galaxies whose present star-formation rates 
are much lower than their past average. A survey of 
star-forming galaxies at z ~ 1 (such as ours) will be only 
negligibly affected by their omission. Of more potential 
concern are galaxies at the opposite extreme, galaxies 
with present star-formation rates much higher than their 
past average. The spectra of these galaxies are dom- 
inated by light from massive O stars — by stars whose 
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atmospheres are too hot to contain significant amounts 
of neutral Hydrogen — and consequently they do not have 
Balmer breaks. Since a Balmer break begins to be dis- 
cernible when a star-formation episode has lasted longer 
than the typical ~ 10 7 yr lifetime of an O star, and since 
most star formation in the local universe (e.g., Heckman 
1997) and at z ~ 3 (e.g., Papovich et al. 2001; Shap- 
ley et al. 2001) occurs in episodes that last substantially 
longer than 10 7 yr, we suspected that our reliance on the 
Balmer-break in our selection criteria would not cause 
us to miss significant numbers of star-forming galaxies 
at z ~ 1. 

In order to obtain empirical support for our proposed 
selection criteria (equation 0, we began in August 1995 
to obtain GTZi images in fields with completed or ongo- 
ing magnitude- limited redshift surveys. Four fields were 
chosen: the 00^53 field of Cohen et al. (1999), the Hubble 
Deep Field (Williams et al. 1996; Cohen et al. 2000), the 
14 h 18 field of Lilly et al. (1995), and the 22 h 18 field of 
Cowie et al. (1996) and Lilly et al. (1995). Spectroscopic 
rcdshifts have been published for 1312 objects in these 
fields combined. The published coordinates for these ob- 
jects were sufficient for us to easily and unambiguously 
identify 1138 of them with objects in our images. 738 
of the objects are in the HDF-N, where spectroscopic 
follow-up has been especially deep and complete; their 
colors and redshifts are shown in figure One can see 
that the bulk of known galaxies at z ~ 1 have GTZi colors 
that satisfy our selection criteria. 



Redshifts in HDF-N 




0.5 



0.5 1 1.5 

z 

Fig. 5. — Observed redshifts of objects with different GTZi colors. 
The data are taken from various magnitude-limited surveys in the 
HDF-N; see text. Known galaxies at z ~ 1 have colors similar to 
our expectations from figure El The solid line shows the z ~ 1 
selection criteria of equationlll 

The galaxies at 0.85 < z < 1.15 that do not satisfy 
our GTZi selection criteria can be crudely grouped into 
three classes. Please refer to figure [SJ First, there are 
the handful of galaxies with TZ — i ~ 0.9, G — TZ ^, 2.0. 



These galaxies have rest-frame colors nearly identical to 
those of local ellipticals (see figure 0J. They are objects 
that formed the bulk of their stars in the past and are 
no longer forming stars at a significant rate. Their ab- 
sence from a color-selected survey of star-forming galax- 
ies is expected and harmless. Second, there are objects 
with colors identical to those of low redshift galaxies but 
with reported spectroscopic redshifts z ~ 1.0. Some of 
these objects may have unusual star-formation histories 
or large photometric errors or exceptionally strong emis- 
sion lines, and some may have incorrectly measured spec- 
troscopic redshifts. Third, there are objects with mea- 
sured colors that lie just outside the Balmer-break color 
selection window. These galaxies may have been scat- 
tered out of our selection window by photometric errors, 
which are typically ~ 0.2 magnitudes in both G — 1Z and 
7Z — i. Their absence from a Balmer-break selected sur- 
vey can be largely corrected with statistical techniques 
described in Adelberger (2002,2004). 

Figure[f)|presents the data of figure|3]in a way that may 
be easier to grasp. We produced one alternate realization 
of figurc[5]s data by adding to each galaxy's GTZi colors a 
Gaussian deviate with standard deviation a = 0.1 that is 
similar to the color's measurement uncertainty. After re- 
peating this procedure numerous times, we concatenated 
the alternate realizations into a large list oiG — TZ, 72. — i, 
z triplets, and then calculated, for each point in the GTZi 
plane, the fraction of galaxies with those colors in the 
large list that had redshift 0.85 < z < 1.15. Figure HO 
shows the result. As anticipated, the portion of the GTZi 
plane that is dominated by star- forming galaxies at z ~ 1 
is approximately described by equation^ One exception 
is the region near G — TZ^ 2.8, TZ — i^ 0.9, which is pop- 
ulated primarily by galaxies at z ~ 1 but lies outside our 
selection window. It would be easy to modify the window 
to include this region, but, as mentioned above, galaxies 
at z <~ 1 with these colors do not account for much of the 
star-formation density. We chose to ignore them. Read- 
ers interested in stellar mass rather than star-formation 
rate might choose differently. 

Ideal color-selection criteria would be perfectly com- 
plete and perfectly efficient; they would be satisfied by 
every galaxy in the targeted redshift interval and only by 
galaxies in the targeted redshift interval. Figure [5] shows 
that in practice the goals of completeness and efficiency 
are incompatible. Photometric errors and intrinsic vari- 
ations in the spectra of galaxies cause galaxies inside the 
redshift interval 0.85 < z < 1.15 to have a wide range 
of GTZi colors. In some cases these colors are identical 
to those of galaxies at other redshifts. If we wanted our 
color-selected survey to be as complete as possible, we 
would want to make our selection box very large so that 
it would include even galaxies with large photometric er- 
rors or abnormal spectral shapes, but this improvement 
in completeness would come at the price of admitting 
more galaxies at the wrong rcdshifts and it would there- 
fore decrease our efficiency. To quantify how closely our 
selection criteria satisfied the conflicting goals of com- 
pleteness and efficiency, we used the magnitude-limited 
surveys discussed above to calculate two quantities, a, 
the completeness coefficient, is equal to the fraction of 
galaxies at 0.85 < z < 1.15 in the magnitude limited sur- 
veys whose colors satisfied equation j3, the efficiency 
coefficient, is equal to the fraction of magnitude- limited 



6 



Adelberger et al. 



Bel. Num. Density, all sources Fraction with 0B5<z<1.15 Rel. Num. Density, 0.8S<2<1.15 
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Fig. 6. — Left panel: The distribution of colors for all objects 
in deep GlZi images of the HDF-N. Middle panel: The proba- 
bility that an object has 0.85 < z < 1.15 as a function of GIZi 
color. Right panel: The distribution of GTZi colors for objects with 
0.85 < z < 1.15 in the HDF-N. Superimposed are the colors of 
the 16 ISO 15/^m sources with 0.84 < z < 1.15 in this field (see 
Adelberger & Steidel 2000 and references therein; 41 of the 46 ISO 
sources in this field have measured spectroscopic redshifts, so it is 
unlikely that there is a large number of spectroscopically uniden- 
tified ISO sources at z ~ 1 with colors significantly different from 
these). These extremely dusty sources are among the most rapidly 
star-forming galaxies at z ~ 1, but their optical colors are suffi- 
ciently similar to those of more normal galaxies that they satisfy 
our selection criteria (equation solid line in all panels). There is 
little evidence that our UV-based selection criteria bias our samples 
against the dustiest galaxies, a point that has been more labori- 
ously made by Adelberger & Steidel (2000) and Adelberger (2001). 



In the fall of 1997 we began to add occasional z ~ 1 
galaxy candidates to our Lyman-break galaxy slitmasks. 
By the spring of 1999 we had settled on our selection 
criteria (equation ^) and were devoting entire slitmasks 
to Balmer-break galaxies. Figure shows the overall 
redshift distribution of the objects observed to date. Ex- 
cluding the handful of galaxies with z < 0.4, whose un- 
expected colors usually resulted from abnormally strong 
nebular emission lines, the mean redshift is (z) = 1.02 
and the r.m.s. is o~ z = 0.15. 
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Fig. 7. — The redshift histogram the spectroscopically observed 
sources in our fields whose colors satisfied the selection criteria of 
equation 111 



survey objects in the selection box of equation whose 
measured redshift satisfied 0.85 < z < 1.15. 

Table [21 lists a and /?, with and without weighting by 
the galaxies' apparent luminosities through various fil- 
ters. Our completeness depends on wavelength. Samples 
selected through equation ^ ar e especially complete for 
the bluest galaxies: the /7 n -weighted column shows that 
approximately 85% of the U n (rest-frame 1800A) lumi- 
nosity density detected at 0.85 < z < 1.15 in magnitude- 
limited surveys is produced by galaxies that satisfy the 
selection criteria. The completeness falls towards redder 
wavelengths, where the total z ~ 1 luminosity density 
receives larger contributions from older galaxies whose 
spectra are less dominated by star formation, but even at 
rest-frame 4500A (observed z-band) approximately 70% 
of the luminosity density at 0.85 < z < 1.15 is produced 
by galaxies whose colors satisfy equation ^ Because our 
color-selected catalogs extend to magnitudes significantly 
fainter than those of the magnitude-limited surveys, the 
detected luminosity density at z ~ 1 in our survey is far 
higher than in the magnitude-limited surveys. The com- 
pleteness fractions above apply only to the magnitude 
range 1Z £ 24 where the surveys overlap. 



4.2. Balmer-break selection at higher redshifts 

In principle it would be easy to use similar color- 
selection criteria to find galaxies at almost arbitrarily 
high redshifts. In practice there is little reason to pur- 
sue this selection strategy beyond the redshift z ~ 1.5 
where the Balmer-break leaves the optical window; near 
this redshift bluer spectral features are beginning to enter 
the optical window and by exploiting these one can con- 
tinue to take advantage of well-developed CCD detector 
technology. 

Little thought is required to extend the Balmer-break 
selection to z ~ 1.5. Consider figure El for example, 
which shows the expected GIZz colors of galaxies at dif- 
ferent redshifts. The color cuts 

K - z > 0.8(G - Tl) + 0.3 G-K< 1.8 (2) 

isolate galaxies at 1.0 < z < 1.5 from the foreground 
and background populations. Figure El calculated in an 
analogous manner to figure El confirms that objects in 
the HDF with these colors tend to lie at redshifts 1.0 < 
z < 1.5. Table [3 lists the completeness parameters a 
and (3 for these selection criteria. Approximately 80% of 
the known luminosity density at redshifts 1.0 < z < 1.5 
and observed wavelengths 3000 < A < 9500A in the HDF 
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Table 2. Color-selection Efficiency, 1 ;$ z ;$ 1.5 



Name 


a, p 


OtU,PU 


Cti,Pi 


CKz, /5z 


Reference 


JV can d 


AT f 


0.85 < z < 1.15: 
















00^53 


0.73,0.55 


0.95,0.55 


0.77,0.55 




Cohen et al. 1999 


34 


168 


HDF-N 


0.67,0.55 


0.81,0.57 


0.67,0.49 


0.65,0.51 


Cohen et al. 2000 


166 


676 


14 h 18 


0.53,0.67 


0.60,0.56 


0.57,0.64 




Lilly et a2. 1995 


15 


123 


22^18 


0.65,0.46 


0.87,0.67 


0.67,0.50 




Cowie et al. 1996 


21 


171 


median 


0.66,0.55 


0.84,0.57 


0.67,0.53 


0.65,0.51 








1.0 < z < 1.5: 
















HDF 


0.80,0.53 


0.84,0.60 


0.78,0.53 


0.77,0.54 


Cohen et al. 2000 


100 


676 



a The fraction of galaxies with 7Z > 21 at the redshift of interest whose colors satisfy our proposed selection criteria («) and the fraction 
of objects with TZ > 21 satisfying our selection criteria whose redshift lies in the desired range (/3). 
b a and /3 recalculated after assigning each galaxy a weight proportional to its apparent U n luminosity. 
c a and f) recalculated after assigning each galaxy a weight proportional to its apparent i luminosity. 
d ct and /3 recalculated after assigning each galaxy a weight proportional to its apparent z luminosity. 
c Number of spectroscopic redshifts in redshift range of interest (TZ > 21) 
f Total number of spectroscopic redshifts (TZ > 21) 



is produced by galaxies whose colors satisfy equation [21 
These numbers should be treated with some caution, 
since the incompleteness of the HDF magnitude-limited 
spectroscopy strongly skews the observed redshift distri- 
bution towards the lower end of the range 1.0 < z < 1.5. 
Nevertheless we hope to have demonstrated that Balmer- 
break selection allows one to create reasonably complete 
catalogs of star-forming galaxies throughout the redshift 
range 1.0 ;$ z £ 1.5. 
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Fig. 8. — The expected GTZz colors of stars and galaxies. The 
symbols are as in figure |I] except here large dots mark the colors 
of model galaxies with E(B - V) = 0.15 at Z = 1.0, 1.1, 1.2, 1.3, 
1.4, 1.5 and the solid line shows the selection criteria of equationl2l 



5. COLOR SELECTION AT 1.9 < Z < 2.7 

Figure shows that the absence of a strong break in 
the optical window is a distinguishing characteristic of 




0.2 0.4 0.6 0.8 1 



Fig. 9. — Similar to figure except results are shown for GTZz 
selection of galaxies with 1.0 < z < 1.5 rather than GTZi selection 
of galaxies with 0.85 < z < 1.15. 



galaxies at z ~ 2. Ruling out the existence of a break re- 
quires photometry through at least the five filters shown 
in the figure, but we wanted to devise selection crite- 
ria that required imaging through only three. We chose 
to use the U n GlZ filters for reasons of convenience; our 
ongoing Lyman-break galaxy survey was producing nu- 
merous deep U n G1Z images that we wanted to use for 
other purposes. 

Since the targeted redshift range 1.9 < z < 2.7 is sim- 
ilar to the redshift range 2.6 < z < 3.4 of the Lyman- 
break survey, we derived our initial estimates of the U G1Z 
colors of galaxies at 1.9 < z < 2.7 from the observed col- 
ors of Lyman-break galaxies. For this we used a sample 
of 70 Lyman-break galaxies that had spectroscopic red- 
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shifts and measured photometry through the U n GlZJK s 
bandpasses. 37 of the galaxies were taken from Shapley 
etal. (2001) and 33 from Papovich et al. (2001). Follow- 
ing the approach outlined in those papers, each galaxy's 
photometry was fit with model spectral energy distribu- 
tions (SEDs) that had a range of star-formation histories 
and dust reddenings. We found the best-fit SED for ev- 
ery galaxy in the sample, redshifted the best-fit SEDs to 
z = 1.5, z = 2.0, and z = 2.5, and calculated their U n GlZ 
colors as described in § [21 This produced an estimate of 
the U n GlZ colors that each Lyman-break galaxy would 
have had if its redshift were z ~ 2 rather than z ~ 3 
(figure ITU)l . Star-forming galaxies at z ~ 2 should have 
U n GlZ colors similar to these. 




I I I I I I I I I I I I I I I I I 

-0.5 0.5 1 

G - % 

Fig. 10. — The expected UnGTZ colors of stars and of star- 
forming galaxies at various redshifts. Circles and squares mark 
the colors that would result from shifting the observed SEDs of 
Lyman-break galaxies at z ~ 3 to z = 1.5, 2.0, 2.5. Data are from 
Papovich et al. (2001; HDF) and Shapley et al. (2001; ground- 
based). Curved tracks show the expected colors of model galaxies 
(§0 a t lower redshifts, starting at z = and increasing clockwise 
to z = 1.5. Crosses mark redshift intervals of dz = 0.1. The Sb 
and Sc tracks move off the plot at 2 ~ 0.2—0.3 and re-enter at 
z ~ 0.7-1.0. Stars represent Galactic stars from Gunn & Stryker 
(1983). The trapezoid is our selection window (equation [3} for 
galaxies with 1.9 z ^ 2.7. 

Useful color-selection criteria must not only find galax- 
ies at the targeted redshifts but also avoid those at other 
redshifts. Galaxies at significantly higher redshifts will 
have extremely red U n — G colors due to the strong Ly- 
man break. They are unlikely to be confused with galax- 
ies at z ~ 2. The tracks on figure show the col- 
ors of model galaxies at lower redshifts. We considered 
the templates Im, Sb, Sc from § [21 and reddened each 
to E(B - V) = 0.15 with dust that followed a Calzetti 
(1997) extinction curve. At 0.3 ^ z ^ 1.0 galaxies have 
red G — 1Z colors due to the redshifted Balmer/4000A 
breaks and are easily distinguished from galaxies at 
z ~ 2. At lower and higher redshifts the potential for 



confusion is greater. Young galaxies (type Im) at z < 0.3 
and 1.5 < z < 1.9 have U n GlZ colors sufficiently simi- 
lar to those of galaxies at 1.9 < z < 2.7 that a clean 
separation is impossible given U n GlZ images alone. Our 
selection criteria were designed in large part to mitigate 
this contamination. 

To learn how to design a selection window that would 
exclude as many interloping galaxies as possible, we be- 
gan in the fall of 1996 to obtain spectra of objects with 
colors similar to those we expected for galaxies at red- 
shift z ~ 2. Figure ITTI shows how galaxies' U n G1Z colors 
depend on their redshifts. In addition to the galaxies 
whose redshifts we measured for this purpose, the fig- 
ure includes all galaxies from the Lyman-break survey of 
Steidel et al. (2003) and all galaxies from the magnitude- 
limited surveys and our Balmer-break survey described 
in §0] FigureslT^landrLUpresent these data in a way that 
may be easier to understand. Both are constructed in a 
similar manner to figure |H1 Figure IT2"1 shows the typical 
U n GlZ colors of galaxies at z < 1. These were the galax- 
ies that we hoped to exclude from our sample. Figure IT51 
shows as a function of U n GlZ color the relative num- 
ber density of sources (left panel), the fraction of objects 
that lie in the redshift range of interest 1.9 < z < 2.7 
(middle panel), and the implied relative number density 
of sources in the same redshift range (right panel). Be- 
cause no existing magnitude-limited surveys contain sig- 
nificant numbers of galaxies at these redshifts, the right 
panel was calculated by multiplying the left and center 
panels together. 



Spectroscopically Observed Galaxies 




0.5 1 1.5 




z 



Fig. 11. — Spectroscopic redshifts of objects with different 
U n G1Z colors. The hooked shaped is produced by the transition 
from Balmer/4000A absorption in the spectra of galaxies at z ^ 1 
to Lyman absorption in the spectra of galaxies at z 2.5. Galax- 
ies at intermediate redshifts have neutral U n G1Z colors; those with 
2.0 ^ z ^ 2.5 tend to lie within the region enclosed by the trape- 
zoid (equation 13*1. 

The following selection criteria were inspired by the 
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probability of satisfying the criteria, and by our desire to 
leave no gap between these criteria and the Lyman-break 
selection criteria of Steidel et al. (2003): 



G 



ft > -0.2 

G>G-K 



0.2 



[BX] 



G-TZ<0.2{U n -G)- 
U n -G<G -TZ +1.0. 



0.4 



(3) 



These criteria are called "BX" in our internal naming 
convention and may be referred to by that name in sub- 
sequent publications. 

Figure ITU shows the redshift distribution of randomly 
selected objects whose colors satisfy these criteria. Ex- 
cluding sources with z < 1, the mean redshift is z = 2.22 
and the standard deviation is a z = 0.34. 



0.2 0.4 0.6 0.8 
fraction with z<l 



Fig. 12. — The fraction of sources that have spectroscopic red- 
shift z < 1 as a function of U n G1i color. The solid lines enclose 
the Lyman-break selection window of Steidel et al. (2003) and the 
two 1.4 < z < 2.7 selection windows discussed below (equations 151 
and El. 
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number density 
23.5<«<25.5 



Implied relative 
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of galaxies with 
1.9<z<2.7 
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Fig. 13. — Left panel: The distribution of colors for all objects 
with K > 23.5 in deep U n GTZ images of the HDF-N. Middle panel: 
The observed probability that an object with TZ > 23.5 has redshift 
1.90 < z < 2.70 as a function of U„G7Z color. Right panel: The 
product of the left and middle panels. This shows the implied 
colors of objects with 1.90 < z < 2.70. The trapezoid in all panels 
shows the selection criteria of equation 151 



shape of the contours on these plots, by our requirement 
that galaxies of all LBG-like spectral types have a finite 



200 



150 - 



^100 




Fig. 14. — The observed redshift distribution of objects whose 
colors satisfy equation [3] Stars are included in the bin containing 
Z = 0; they make up 4.6% of the sample but only 0.9% of the 
sample fainter than TZ = 23.5. 



As could be anticipated from figure the Balmer 
break in low redshift galaxies can sometimes be confused 
with the Lyman-a forest absorption in the spectra of 
galaxies at z ~ 2. The resulting contamination of our 
sample is severe at magnitudes TZ < 23 but negligible by 
TZ = 25.5. Restricting the sample to TZ > 23.5 provides 
a crude but effective way of eliminating low-redshift in- 
terlopers. 

One can roughly estimate the completeness coefficients 
a and j3 for the sample as follows. Let P(x,y) be the 
probability that a randomly chosen galaxy with 23.5 < 
TZ < 25.5 and colors G — TZ = x and U n — G = y has a 
redshift that lies in the range 1.9 < z < 2.7, let n(x, y) be 
proportional to the observed number density of galaxies 
with 23.5 < TZ < 25.5 that have colors G — TZ = x and 
U n — G = y, and let I(x, y) be equal to 1 if the color G — 
TZ — x, U n — G = y satisfies equation[3]and to otherwise. 
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Then the probability than an object with 23.5 < ft < 
25.5 at 1.9 < z < 2.7 will satisfy our selection criteria is 



dx dy P(x, y)n(x, y)I(x, y) 




dxdy P(x,y)n(x,y) 
(4) 



and the probability that an object with 23.5 < ft < 25.5 
that satisfies our selection criteria will lie at 1.9 < z < 2.7 
is 




0= dxdy P(x,y)n(x,y)I(x,y) 



(5) 

The functions P(x, y) and n(x, y) are shown in figure IT31 
Numerically integrating equations 0] and [S] yields the es- 
timates a — 0.64, — 0.70. Roughly two-thirds of galax- 
ies with 1.9 < z < 2.7 satisfy our selection criteria, and 
roughly two-thirds of the galaxies that satisfy our selec- 
tion criteria lie at 1.9 < z < 2.7. The result could have 
been anticipated to a large extent from the redshift his- 
togram shown in figure ITU 

6. COLOR SELECTION AT 1.4 < Z < 2.1 

Our approach towards defining selection criteria at 
these redshifts differs little from our approach at higher 
redshift. As redshift decreases from z = 2, galaxies be- 
come bluer in U n — G because a smaller fraction of their 
observed-frame U n emission is absorbed by the Lyman- 
a forest. Otherwise their optical colors are largely un- 
changed. One would expect galaxies with 1.4 < z < 2.1 
to lie just below our 1.9 < z < 2.7 selection window 
(equation^ in the U n GlZ plane. We began exploratory 
spectroscopy of galaxies in this part of the U n GlZ plane 
in the fall of 1997. Observations continued sporadically 
until the spring of 2003. Figure [TBI shows that these ob- 
servations largely conformed to our expectations. The 
following selection criteria were inspired by the shape of 
the contours on this plot, by our requirement that galax- 
ies of all LBG-like spectral types have a finite probability 
of satisfying the criteria (see figurc lTU|) , and by our desire 
to leave no gap between these criteria and the selection 
criteria of equation 



G-ft>-0.2 
U n -G>G -ft -0.1 

G-ft<0.2([/ n -G)+0.4 
U n -G<G -ft + 0.2. 



[BM] 



(6) 

These criteria are called "BM" in our internal naming 
convention and may be referred to by that name in sub- 
sequent publications. 

The completeness coefficients a and /? for these selec- 
tion criteria, estimated with the approach of equations 2] 
and El are listed in tableau Roughly 42% of galaxies with 
1.4 < z < 2.1 and ft < 25.5 satisfy the selection criteria, 
and roughly 46% of the objects with ft < 25.5 that sat- 
isfy the criteria are galaxies with 1.4 < z < 2.1. FigurelTfil 
shows the observed redshift distribution of the sources 
whose colors satisfied equation|S| Excluding sources with 
z < 1, the mean redshift is z = 1.70 and the standard de- 
viation is <t z = 0.34. As figure El shows, galaxies' U n GlZ 
colors do not change by a large amount between z = 1.0 
and z = 1.5, and as a result the redshift distribution has 
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dxdy n(x,y)I(x,y). 



Fraction with 
1.5<z<2.0 



Implied relative 
number density 
of galaxies with 
1.4<z<2.1 
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Fig. 15. — Similar to figurelBI except the targeted redshift range 
is 1.4 < z < 2.1, the TZ > 23.5 criteria has been removed, and the 
trapezoids show the selection criteria of equation IB1 



Table 3. Color-selection Efficiency, 1.4 ^ 2 ^ 2.7 



Redshift 


o,/3 a 


1.4 < z < 2.1: 


0.42,0.46 


1.9 < z < 2.7: 


0.64,0.70 



a The estimated fraction of galaxies at the redshift of interest 
whose colors satisfy our proposed selection criteria (a) and the 
estimated fraction of objects satisfying our selection criteria whose 
redshift lies in the desired range (/3). 



a significant tail extending to z ~ 1. This low-redshift 
tail can be eliminated to a large extent, if z band pho- 
tometry is available, by excluding from the sample any 
objects whose colors satisfy the 1.0 < z < 1.5 selection 
criteria of equation |3 See figure 1171 Subsequent papers 
may refer to this combination of color-selection criteria 
as "BMZ." 



7. 



U N GTZ color selection at any redshift 

1< Z < 3 



The previous two sections presented color-selection cri- 
teria tuned to the two arbitrary redshift ranges 1.4 < z < 
2.1 and 1.9 < z < 2.7. Some projects may require a sam- 
ple of galaxies with a similar but slightly different range 
of redshifts, e.g., 1.7 < z < 2.3. Minor adjustments to 
the selection criteria we have presented can tune them 
to this redshift range or others. To help readers esti- 
mate how to adjust our criteria to produce samples with 
a desired range of redshifts, we show in figure El the 
observed median redshift as a function of t/„Gft color 
of galaxies with ft > 23.5 in our spectroscopic sample. 
Roughly aligning the edges of a selection box with this 
plot's contours will produce reasonably good selection 
criteria tuned to arbitrary redshifts within the interval 
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Fig. 16. — The redshift distribution of randomly selected objects 
whose colors satisfy equation IB1 
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Fig. 17. — Light histogram: the redshift distribution of spec- 
troscopically observed objects in the HDF-N whose colors satisfy 
equation [H] The stellar contamination rate is about 1%. Dark his- 
togram: the redshift distribution of HDF-N objects whose colors 
satisfy equation HJ but not equation 121 



1 ^ z ,$ 3. Some spectroscopic follow-up will be required 
to verify the median redshift of the sources and to place 
limits on the sample's contamination by stars and low- 
redshift galaxies. 

8. SUMMARY AND DISCUSSION 



O 

C 



0.5 1 

G-31 



1 1 1 1 L i 1^1 1 J 1 


1.5 2 2 


5 3 3 


r i i i i i i i i 
ii 


r i t i i i 
i i i i i 



0.2 



0.4 0.6 
median z 



0.8 



Fig. 18. — Median redshift of spectroscopic sources as a function 
of U n GTZ color. Galaxies from the U n G1Z-selected spectroscopic 
sample described in this paper, from the Lyman-break survey of 
Steidel et al. (2003), and from the magnitude-limited surveys listed 
in §g]were placed into bins with A(U„ — G) = 0.1, A(G-TZ) = 0.1 
and the median redshift was calculated for each bin. The resulting 
array of median redshift vs. U n G1Z color was then smoothed by a 
Gaussian with ajj n —c = CG-TZ = 0-1 to produce this plot. The 
slow change of median redshift with color makes it easy to devise 
selection criteria that target arbitrary redshift ranges within the 
larger span 1 < z < 3. 



It is often asserted that studying galaxies at redshifts 
1 < z < 3 will be tremendously difficult. The opti- 
cal spectra of galaxies near the middle of this redshift 
range contain neither the strong spectral breaks that are 
sometimes thought to be necessary for effective photo- 
metric selection nor the strong emission lines Ly-a or 
[OIIJA3727 that are sometimes thought to be crucial for 
spectroscopic identification. The supposed difficulty of 
galaxy observations causes many to refer to this redshift 
range as the spectroscopic desert or high place of sacrifice 
(Bullock et al. 2001). 

We are sceptical. The optical colors of galaxies at red- 
shifts 1 < z < 3 are distinctive. We showed in sec- 
tions 2)13 and H3 that these galaxies are easy to locate in 
deep images using simple color-selection criteria (equa- 
tion for 0.9 < z < 1.1, equation for 1.0 < z < 1.5, 
equation [3] for 1.9 < z < 2.7, and equation for 
1.4 < z < 2.1). Once they have been found their red- 
shifts are no harder to measure than those of comparably 
bright galaxies at slightly higher or lower redshift. This 
is due in large part to the strength of their interstellar 
absorption lines. With an appropriately chosen spectro- 
scopic set-up, redshifts are as easy to measure from ab- 
sorption lines as from the Lyman-a or [OII]A3727 emis- 
sion lines. The point is illustrated by figure 1191 during 
our surveys we obtained as many redshifts per hour of 
observing time at 1.5 < z < 2.5 as at the higher redshifts 
where Lyman-a is more easily observed. That was pos- 
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sible only because we knew the approximate redshifts of 
our sources in advance and could optimize our choice of 
spectrograph and its configuration accordingly. We could 
not have measured a redshift for many of the galaxies at 
1.5 < z < 2.5 with a spectrograph that lacked the good 
UV throughput of LRIS-B. See Steidel et al. (2004) for 
a more detailed discussion. 




& 10° 



large and ongoing spectroscopic surveys of Davis et al. 
(2003) and Le Fevre et al. (2003). 
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Fig. 19. — Scaled redshift distributions of objects selected with 
equations [5] and [3] and with the z ~ 3 Lyman-break selection 
criteria of Steidel et al. (2003). The actual redshift histogram of 
each sample was divided by the spectroscopic observing time de- 
voted to the sample. All histograms were then multiplied by the 
same arbitrary constant and fit with a cubic spline. The area under 
each curve is proportional to the number of spectroscopic redshifts 
we measured per hour with the LRIS spectrograph. Galaxy spec- 
troscopy is not much harder at 1 < z < 3 than at higher or lower 
redshifts — as long as the spectrograph has good UV throughput. 

Two to three nights spent imaging a single field with a 
UV-sensitive 40' x 40' camera on a 4m telescope is suffi- 
cient to detect ~ 2 x 10 4 galaxies that satisfy one set of 
the 1 < z < 3 color-selection criteria that we have pre- 
sented. A major benefit of color-selected spectroscopy is 
that it lets one draw statistically significant conclusions 
from this large number of high-redshift galaxies without 
having to measure a redshift for every one. For example, 
the ~ 1000 spectroscopic redshifts we have measured for 
objects whose colors satisfy equationHtell us, with high 
precision, what the redshift distribution must be for the 
large ensemble of objects in the 40' x 40' image whose 
colors satisfy equation Q] One can use this knowledge 
to derive the luminosity function or spatial clustering 
strength of galaxies at z ~ 1 from the list of photo- 
metric candidates alone. Because the redshift distribu- 
tion of photometric candidates does not depend strongly 
on magnitude for faint galaxies (figure EPjl ■ it should be 
possible to use purely photometric observations to learn 
about the properties of the numerous high-redshift galax- 
ies that are too faint for spectroscopy. This may be best 
use of our z ~ 1 color-selection criteria, since the bright- 
est galaxies at this redshift are routinely detected in the 



Fig. 20. — Redshift distributions as a function of apparent mag- 
nitude. Shaded regions show the mean ±lcr redshift range for ob- 
jects of different apparent magnitudes whose colors satisfy (from 
bottom to top) equations HI 151 13l and the z ~ 3 Lyman-break selec- 
tion criteria of Steidel et al. (2003). The lc interval was calculated 
after excluding interlopers with z < 0.4 from the z ~ 1.0 sample 
and interlopers with z < 1 from the other samples. The number 
in each box indicates the fraction of sources that were excluded. 
This fraction depends strongly on magnitude for some samples, but 
otherwise the redshift distributions are insensitive apparent magni- 
tude. It is therefore reasonable to assume, when the interloper frac- 
tion is small, that faint photometric candidates will have the same 
redshift distribution as brighter candidates with spectroscopic red- 
shifts. The weak trend of higher redshifts for (apparently) fainter 
objects results from the change in luminosity distance. 

The weakness of color-selected surveys is that not all 
galaxies at the targeted redshifts will satisfy the adopted 
selection criteria. A high level of completeness at the 
targeted redshifts can be obtained only at the price of 
admitting large numbers of galaxies at the wrong red- 
shifts. With the simple selection criteria we presented, 
observers will find more than ~ one-half of the galaxies 
brighter than the magnitude limit at the redshift of in- 
terest and will waste no more than ~ one-half of their 
observing time observing objects at other redshifts. 50% 
completeness is not ideal, but the problem is less severe 
than one might imagine. First, a star-forming galaxy at 
1 < z < 3 that does not satisfy one of our color-selection 
criteria will likely satisfy another. Consider galaxies with 
redshift 1.8 < z < 2.2 and magnitude 1Z < 25.5, for ex- 
ample. Equation 0] implies that only 63% of them will 
satisfy the 1.9 < z < 2.7 selection criteria of equation|3| 
But an additional 26% will satisfy the 1.5 < z < 2.0 se- 
lection criteria of equation El and an additional 5% will 
satisfy the Lyman-break galaxy selection criteria of Stei- 
del et al. (2003). A total of 94% of photometrically de- 
tected galaxies with 23.5 < TZ < 25.5 and 1.8 < z < 2.2 
will satisfy one of the U n GlZ selection criteria we have 
presented. By conducting redshift surveys with each of 
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these criteria, one can reduce the incompleteness to a 
reasonably low level. Second, even if a survey adopts 
only one of our selection criteria, much of the incom- 
pleteness can be corrected in a statistical sense. The in- 
completeness is largely due to the photometric errors in 
our color measurements, which are not small compared 
to the size of our color-selection window. Many galax- 
ies whose true colors lie inside our selection window will 
have measured colors that lie outside the window; many 
galaxies whose measured colors lie inside will have true 
colors that lie outside. The result is a broad, bell-shaped 
redshift histogram rather than a boxcar extending from 
the minimum to the maximum targeted redshift. Be- 
cause photometric uncertainties are easy to characterize 
with Monte-Carlo simulations, their contribution to our 
incompleteness is easy to understand and correct. Adel- 
berger (2002,2004) explains in detail. 

In any case, all observational strategies require some 
compromise between efficiency and completeness, and 
the compromises of color selection do not look bad com- 
pared to the alternatives. Suppose one were interested 
in studying galaxies at 0.85 < z < 1.15. By obtain- 
ing a spectrum of every galaxy brighter than 1Z = 24.0 
in an image one would be certain to produce a statis- 
tically complete sample of galaxies at z ~ 1 to this 
magnitude limit, but over 80% of the observing time 
would have been spent on objects at the wrong redshifts 
(e.g., Cohen et al. 2000). If one were willing to toler- 
ate 15% incompleteness in observed U n luminosity den- 
sity, one could cull spectroscopic targets with the color 
criteria of equation ^ an d reduce the required observ- 
ing time by more than a factor of three (see table |3J) . 
As redshift increases the benefits of color-selection be- 
come more obvious. Only one object out of 50 in the 
magnitude limited survey of Cohen et al. (2000) had a 
redshift 1.9 < z < 2.7. If one wished to study galax- 
ies at these redshifts with a magnitude limited survey 
the required observing time would be 33 times longer 
than if one color-selected targets with equation|2] (see ta- 
bic A project that could be completed in 1-2 years 
with color-selection would require a lifetime of observ- 
ing with standard magnitude- limited techniques. The 
~ 30% incompleteness of color-selected surveys does not 
seem a high price to pay. 6 

6 A magnitude-limited survey conducted with a blue-optimized 
spectrograph would find a higher fraction of sources with 1.9 < z < 
2.7 than Cohen's 1 in 50, so the comparison is somewhat unfair. 
In practice, however, few would chose to undertake a magnitude- 
limited survey with a blue spectrograph because it would make 
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Magnitude-limited optical surveys are not the only al- 
ternative to color-selected optical surveys. Many have 
advocated finding galaxies at 1 < z < 3 with photome- 
try outside of the optical window. We chose to develop 
criteria that relied solely on optical photometry because 
ground-based optical imagers offer an unrivaled combi- 
nation of high sensitivity, high spatial resolution, and 
large fields-of-view. This is an advantage that is difficult 
to overcome. Surveys at other wavelengths have their 
strengths — near-IR observations should provide a more 
complete census of older stars (e.g., Franx et al. 2003, 
Rudnick et al. 2003) and far-IR/sub-mm observations 
should yield more reliable estimates of star-formation 
rates — but the fact remains that an investigator given 
two nights on an older 4m optical telescope can create a 
photometric sample of high-redshift galaxies larger than 
all existing samples at other wavelengths combined. It 
seems likely to us that much of what we will learn about 
galaxies at 1 < z < 5 in the coming decade will come 
from large optical surveys selected with color criteria sim- 
ilar to ours. 

Attentive readers may have recognized that the deriva- 
tion of these criteria did not require much thought. That 
is exactly the point. The value of this paper, if any, lies 
not in our photometric selection criteria themselves but 
in the proof that large samples of galaxies at 1 < z < 3 
can be created with trivial techniques that rely solely on 
ground-based imaging and well-developed CCD technol- 
ogy. The spectroscopic desert is a myth. 
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spectra more difficult to identify at the lower redshifts where most 
sources lie. The ability to take full advantage of optimized spec- 
trographs is a non-negligible benefit of color selection. 
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